S3 



ABSTRACT 



Early-type galaxies in the SDSS. I. The sample 
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A sample of nearly 9000 early-type galaxies, in the redshift range 0.01 < z < 0.3, was selected 
from the Sloan Digital Sky Survey using morphological and spectral criteria. This paper describes 
how the sample was selected, presents examples of images and seeing corrected fits to the observed 
■ surface brightness profiles, describes our method for estimating K-corrections, and shows that the 

SDSS spectra are of sufficiently high quality to measure velocity dispersions accurately. It also 
provides catalogs of the measured photometric and spectroscopic parameters. In related papers, 
these data are used to study how early-type galaxy observables, including luminosity, effective 
£T) ' radius, surface brightness, color, and velocity dispersion, are correlated with one another. 
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Subject headings: galaxies: elliptical — galaxies: evolution — galaxies: fundamental parameters 
— galaxies: photometry — galaxies: stellar content 



1. Introduction 

Galaxies have a wide range of luminosities, colors, masses, sizes, surface brightnesses, morphologies, 
star formation histories and environments. This heterogeneity is not surprising, given the variety of physical 
processes which likely influence their formation and evolution, including gravitational collapse, hydrody- 
namics, turbulence, magnetic fields, black-hole formation and accretion, nuclear activity, tidal and merger 
interactions, and evolving and inhomogeneous cosmic radiation fields. 

What is surprising is that populations of galaxies show several very precise relationships among their 
measured properties. The properties we use to describe galaxies span a large "configuration space", but 
galaxies do not fill it. Galaxy spectral energy distributions, when scaled to a fixed broad-band luminosity, 
appear to occupy a thin, one-dimensional locus in color space or spectrum space (e.g., Connolly & Szalay 
1999). Spiral galaxies show a good correlation between rotation velocity and luminosity (e.g., Tully & Fisher 
1977; Giovanelli et al. 1997). Galaxy morphology is strongly correlated with broad-band colors, strengths 
of spectral features, and inferred star-formation histories (e.g., Roberts & Haynes 1994). 

Among all galaxy families, early-type (elliptical and SO) galaxies show the most precise regularities 
(Djorgovski & Davis 1987; Burstcin, Bender, Faber & Nolthenius 1997). Early- type galaxy surface-brightness 
distributions follow a very simple, universal "de Vaucouleurs" profile (de Vaucouleurs 1948). Their spectral 
energy distributions appear to be virtually universal, showing very little variation with mass, environment, or 
cosmic time (e.g., van Dokkum & Franx 1996; Pahre 1998). What variations they do show are measurable and 
precise. Early-type galaxy colors, luminosities, half-light radii, velocity dispersions, and surface brightnesses 
are all correlated (Baum 1959; Fish 1964; Faber & Jackson 1976; Kormendy 1977; Bingelli, Sandage & 
Tarenghi 1984); they can be combined into a two-dimensional "Fundamental Plane" with very little scatter 
(e.g., Dressier et al. 1987; Djorgovski & Davis 1987; Faber ct al. 1987). 

The homogeneity of the early-type galaxy population is difficult to understand if early-type galaxies are 
assembled at late times by stochastic mergers of less-massive galaxies of, presumably, different ages, star 
formation histories, and gas contents, as many models postulate (e.g., Larson 1975; White & Rees 1978; van 
Albada 1982; Kauffmann 1996; Kauffmann & Chariot 1998). It is possible that the homogeneity of early-type 
galaxies points to early formation (e.g., Worthey 1994; Bressan et al. 1994; Vazdekis et al. 1996; Tantalo, 
Chiosi & Bressan 1998); certainly their stellar populations appear old (e.g., Bernardi et al. 1998; Collcss 
ct al. 1999; Trager et al. 2000a, b; Kuntschner et al. 2001). Alternatively, the observable properties of the 
stellar content of early-type galaxies are fixed entirely by the properties of the collisionlcss, self-gravitating, 
dark-matter haloes in which we believe such galaxies lie (e.g., Hernquist 1990). These halos, almost by 
definition, are not subject to the vagaries of gas dynamics, star formation, and magnetic fields; they are 
influenced only by gravity. 

It is essentially a stated goal of the Sloan Digital Sky Survey (SDSS; York et al. 2000; Stoughton et al. 
2002) to revolutionize the study of galaxies. The SDSS is imaging tt steradians of the sky (Northern Galactic 
Cap) in five bands and taking spectra of ~ 10 6 galaxies and ~ 10 5 QSOs. Among the 10 6 SDSS spectra 
there will be roughly 2 x 10 5 spectra taken of early-type galaxies; in fact 10 5 of the spectroscopic fibers are 
being used to assemble a sample of luminous early-type galaxies with a larger mean rcdshift than the main 
SDSS sample (Eisenstein et al. 2001). The high quality of the SDSS 5-band CCD imaging (Gunn et al. 
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1998; Lupton et al. 2001) allows secure identification of early-type galaxies and precise measurements of their 
photometric properties; most spectroscopic targets in the SDSS are detected in the imaging at signal-to-noise 
ratios (S/N) > 100. 

Early-type galaxy studies in the past, for technical reasons, have concentrated on galaxies in clusters at 
low (e.g., J0rgensen, Franx & Kjaergaard 1996; Ellis et al. 1997; Pahre, Djorgovski & de Carvalho 1998a, b; 
Scodeggio et al. 1998; Colless et al. 2001; Saglia et al. 2001; Kuntschner et al. 2001; Bernardi et al. 2002a,b) 
and intermediate rcdshifts (e.g., van Dokkum ct al. 1998, 2001; Kelson et al. 2000; Ziegler et al. 2001). 
Only the large area 'Seven Samurai' (e.g., Faber et al. 1989) and ENEAR surveys (e.g., da Costa et al. 
2000) of nearby early-types, recent work with galaxies in the SBF survey (Blakeslee et al. 2001), and some 
studies at intermediate redshifts by Schade et al. (1999), Treu et al. (1999, 2001a, b) and van Dokkum et 
al. (2001), are not restricted to cluster environments. In constrast, the SDSS is surveying a huge volume 
of the local Universe, so the sample includes early-type galaxies in every environment from voids to groups 
to rich clusters. As of writing, when only a small fraction of the planned SDSS imaging and spectroscopy 
has been taken, the number (~ 9000) of early-type galaxies with well-measured velocity dispersions and 
surface-brightness profiles in the SDSS greatly exceeds the total number in the entire astronomical literature 
to date. 

This is the first of four papers in which we use the SDSS sample to measure the Fundamental Plane 
and other early-type galaxy correlations in multiple bands. This first paper describes how the sample was 
selected and presents the data. Section 2 describes the main properties of the SDSS database. The criteria 
used to select the early- type galaxy sample are described in Section 3. This section also presents a selection 
of images from a range of redshifts, shows de Vaucouleur profile fits to the observed surface brightness 
profiles, and argues that the SDSS photometric pipeline treatment of seeing and sky-subtraction have not 
strongly compromised the estimates of the best-fitting model parameters. It also displays examples of the 
spectra from which we estimate velocity dispersions. Section 4 presents the final catalog of photometric 
and spectroscopic parameters which we use for our subsequent analyses of early- type galaxy properties, and 
discusses some properties of the sample. Many details are relegated to Appendices. Appendix A contains a 
discussion of the various K-corrections we have tried. The way we estimate velocity dispersions is presented 
in Appendix B. A novel method for estimating aperture corrections to the velocity dispersions is discussed 
in Appendix C. Finally, the covariance matrix of the errors is discussed in Appendix D. 

Paper II (Bernardi et al. 2003a) studies correlations between various pairs of observables, such as 
the Faber- Jackson and Kormendy relations. It also presents the luminosity function and its evolution. The 
Fundamental Plane and its dependence on waveband, color, redshift, and environment is studied in Paper III 
(Bernardi et al. 2003b). The co-added spectra of these galaxies are studied in Paper IV (Bernardi et al. 
2003c) . One of the results of that paper is a library of co-added spectra which contains spectra that represent 
a wide range of early- type galaxies. This library is available electronically. These spectra indicate that the 
chemical composition of the early- type galaxy population evolves with redshift. The chemical abundances 
and evolution are then combined with stellar population models to estimate the ages and metallicities of the 
galaxies in our sample. This paper also analyzes correlations with color (e.g., color-magnitude and color-cr 
relations) and discusses the effects of color gradients on measurements of the strength of the correlation 
between color and magnitude. 

Except where stated otherwise, we write the Hubble constant as Ho = 100 h kms -1 Mpc -1 , and we 
perform our analysis in a cosmological world model with (^m,^a,^) = (0.3,0.7,0.7), where and f^A 
are the present-day scaled densities of matter and cosmological constant. In such a model, the age of 
the Universe at the present time is t n = 9.43/i _1 Gyr. For comparison, an Einstein-de Sitter model has 
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(fiM, ^a) = (1,0) and t = 6.52/i _1 Gyr. We frequently use the notation hjo as a reminder that we have set 
h = 0.7. Also, we will frequently be interested in the logarithms of physical quantities. Our convention is to 
set R = log 10 R and V = log 10 a, where R and a are effective radii in h^Q kpc and velocity dispersions in 
km s _1 , respectively. 

2. The SDSS database 

The SDSS project is described in Stoughton et al. (2002). The data we analyze in this paper were 
selected from the SDSS database in summer 2001. At that time the SDSS had imaged ~ 1, 500 square degrees; 
~ 65, 000 galaxies and ~ 8000 QSOs had both photometric and spectroscopic information. The photometric 
and spectroscopic data were taken with the 2.5-m SDSS Telescope at the Apache Point Observatory (New 
Mexico) between 1999 March and 2000 October. Details of the photometric and spectroscopic observations 
and data reduction procedure will be presented elsewhere. Here we briefly summarize. 

2.1. SDSS imaging data 

Images are obtained by drift scanning with a mosaic CCD camera (Gunn et al. 1998) which gives a field 
of view of 3 x 3 deg 2 , with a spatial scale of 0.4 arcsecpix -1 in five bandpasses (u, <?, r, i, z) with central 
wavelengths (3560, 4680, 6180, 7500, 8870A) (Fukugita et al. 1996). The errors in u band measurements are 
larger than the others, so we will only present results in the other four bands. In addition, the photometric 
solutions we use in this paper are preliminary (for details, see discussion of the Early Data Release in 
Stoughton et al. 2002); we use r* rather than r, and similarly for the other bands, to denote this. 

The effective integration time is 54 sec. The raw CCD images are bias-subtracted, flat-fielded and 
background-subtracted. Pixels contaminated by the light of cosmic rays and bad columns are masked. 
Astronomical sources are detected and overlapping sources are de-blended. The data are flux-calibrated by 
comparison with a set of overlapping standard-star fields calibrated with a 0.5-m "Photometric Telescope" 
(Hogg et al. 2001; Smith et al. 2002). The Photometric Telescope is also used for measuring the atmospheric 
extinction coefficients in the five bands. The median effective seeing (the median FWHM of the stellar 
profiles) for the observations used here is 1.5 arcsec. All of this image processing is performed with software 
specially designed for reducing SDSS data (Lupton et al. 2001). The uncertainty in the r*-band zero-point 
calibration is < 0.01 mag; the uncertainty in the sky background subtraction is less than about 1%. 

Surface photometry measurements are obtained by fitting a set of two-dimensional models to the images. 
The model fits account for the effects of seeing, atmospheric extinction, and Galactic extinction (this last 
uses the results of Schlegel, Finkbeiner, & Davis 1998). The SDSS model for seeing is described in Stoughton 
et al. (2002). Briefly, the PSF is expected to have a Fourier Transform oc cxp[— (/crpsF/2) 5//3 ] characteristic 
of Kolmogorov turbulence (e.g., Saglia et al. 1993). This PSF is approximated as the sum of Gaussians (this 
parametrization allows for a substantial reduction in the processing time required to deconvolve the effects 
of seeing) . 

The SDSS photometric pipeline fits two models to the two-dimensional image of each object in each band: 
a pure de Vaucouleurs profile and a pure exponential profile (e.g. Stoughton et al. 2002). Briefly, the photon 
counts are binned into a number of radial and twelve angular bins. Then, two-dimensional de Vaucouleurs 
and exponential models, convolved with the seeing, are fitted to the cumulative binned counts. The pipeline 
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does not include more complicated models (e.g., bulge plus disk) because the computational expense to 
compute them is not justified for the majority of the detected objects. Algorithms which fit bulge and disk 
components (which are desirable to fit on large galaxies) to the surface brightness profiles are currently being 
developed within the collaboration, but are not yet available. We present examples of r*-band images and 
mean surface brightness profiles in Section 3.2.2. 

The SDSS image processing software provides several global photometric parameters, for each object, 
which are obtained independently in each of the five bands. Because we are interested in early-type galaxies, 
we use primarily the following: 

• The ratio b/a of the lengths of the minor and major axes of the observed surface brightness profile. 

• The effective radius (or half-light radius) rd cv along the major axis and 

• the total magnitude mdev; these are computed by fitting a two-dimensional version of de Vaucouleurs 
(1948) r 1 / 4 model to the observed surface brightness profile. (The fitting procedure accounts for the 
effects of seeing — we test the accuracy of this procedure below.) 

• Likelihood parameters deV_L and exp_L that indicate how well the de Vaucouleurs or exponential 
models, when convolved with a model for the seeing, fit the observed light profile. 

• The model magnitude m m ; this is the total magnitude calculated by using the (de Vaucouleurs or 
exponential) model which fits the galaxy profile best in the r*-band. The model magnitudes in the 
other four bands are computed using that r* fit as filter; in effect, this measures the colors of a galaxy 
through the same aperture. 

• The Petrosian magnitude m p is also computed; this is the flux within 2r p , where r p is defined as the 
angular radius at which the ratio of the local surface brightness at r to the mean surface brightness 
within a radius r is 0.2 (Petrosian 1976). 

• The Petrosian radii r 50 and r go ; these are the angular radii containing 50% and 90% of the Petrosian 
light, respectively. 

Although our analyses on early-type galaxy properties presented here and in the companion papers have 
been performed with both the de Vaucouleurs fit parameters and the Petrosian quantities, in most of our 
analyses only the results of the de Vaucouleurs fits are presented. This is because the de Vaucouleurs model 
appears to be a very good fit to the early-type galaxy surface-brightness profiles in the SDSS sample and 
because it is conventional, in the literature on early-type galaxies, to use these quantities. On the other hand, 
for reasons given in Stoughton et al. (2002), unless stated otherwise, galaxy colors are always computed 
using model magnitudes. 

2.2. SDSS spectroscopic data 

The SDSS takes spectra only for a target subsamplc of objects. Target selection criteria are described 
in Stoughton et al. (2002) and Strauss et al. (2002). Spectra are obtained using a multi-object spectrograph 
which observes 640 objects at once. Each spectroscopic plug plate, 1.5 degrees in radius, has 640 fibers, each 
3 arcsec in diameter. Two fibers cannot be closer than 55 arcsec due to the physical size of the fiber plug. 



- 6 - 



Typically ~ 500 fibers per plate are used for galaxies, <~ 90 for QSOs, and the remaining for sky spectra and 
spectrophotometric standard stars. 

Each plate typically has three to five spectroscopic exposures of fifteen minutes, depending on the 
observing conditions (weather, moon); a minimum of three exposures is taken to ensure adequate cosmic ray 
rejection. For galaxies at z < 0.3 the median spectrum S/N per pixel is 16 (see Figure 18 in Appendix B). The 
wavelength range of each spectrum is 3900 — 9000 A. The instrumental dispersion is log 10 A = 10 _4 dex/pixel 
which corresponds to 69 km s _1 per pixel. (There is actually some variation in this instrumental dispersion 
with wavelength, which we account for; see Figure 16 and associated discussion in Appendix B.) The 
instrumental resolution of galaxy spectra, measured from the autocorrelation of stellar template spectra, 
ranges from 85 to 105 km s _1 , with a median value of 92 km s^ 1 . 

A highly automated software package has been designed for reducing SDSS spectral data. The raw data 
are bias-subtracted, flat-fielded, wavelength calibrated, sky-lines removed, co-added, cleaned from residual 
glitches (cosmic rays, bad pixels), and flux calibrated. The spectro-software classifies objects by spectral 
type and determines emission and absorption redshifts. (Rcdshifts are corrected to the heliocentric reference 
frame.) The rcdshift success rate for objects targeted as galaxies is > 99% and errors in the measured redshift 
are less than about 10~ 4 . Once the redshift as been determined the following quantites are computed: 
Absorption-line strengths (Brodie & Hance 1986; Diaz, Terlevich & Terlevich 1989; Trager et al. 1998), 
equivalent widths of the emission lines, and eigen-coefficients and classification numbers of a PCA analysis 
(Connolly & Szalay 1999). Some information about the reliability of the redshift and the quality of the 
spectrum is also provided. We present examples of the spectra in Section 3.3.2. 

The SDSS pipeline does not provide an estimate of the line-of-sight velocity dispersion, er, within a 
galaxy, so we compute it separately for the early-type galaxy sample (see Section 3.3.1). 

3. The sample 

3.1. Selection criteria 

The main goal of this series of papers is to study the properties of early-type galaxies using the main 
galaxy sample of the SDSS database. Therefore, one of the crucial steps in our study is the separation of 
galaxies into early and late types. We want to select objects whose spectra are good enough to compute the 
central velocity dispersion. In addition, because we wish to study the colors of the galaxies in our sample, 
we must not use color information to select the sample. To reach our goal, we have selected galaxies which 
satisfy the following criteria: 

• concentration index rgo/r^o > 2.5 in «*; 

• the likelihood of the de Vaucouleurs model is at least 1.03 times the likelihood of the exponential model; 

• spectra with PCA classification numbers a < —0.1, typical of early-type galaxy spectra (Connolly & 
Szalay 1999); 

• spectra without masked regions (the SDSS spectroscopic pipeline outputs a warning flag for spectra of 
low quality; we only chose spectra for which this flag was set to zero) ; 

• S/N > 10; 
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• redshift < 0.3. 

Section 3.3.1 describes how we estimate the velocity dispersion a for galaxies which satisfy the above 
criteria. We consider velocity dispersion estimates smaller than about 70 km s _1 to be unreliable (see 
Appendix B). Therefore, in the final sample presented in Section 4, galaxies with a < 70 km s _1 have been 
excluded. The results of this and the companion papers are not significantly different if we change the cut-off 
on velocity dispersions to 100 km s _1 . 

About 9,000 objects satisfied all the above criteria. 

As stated earlier, the SDSS pipeline does not output disk-to-bulge ratios from fits to the light profiles. 
The first two requirements above attempt to select profile shapes which are likely to be those of spheroidal 
systems. The i* band measurements tend to be less noisy than g* or z* , so we chose to use the i* band 
estimate of how centrally concentrated the light is. The second requirement reduces approximately to 
requiring that deV_L/ (deV_L + exp_L) > 0.5; that is, the surface brightness profile should be better fit by the 
de Vaucouleurs model than by an exponential. 

The spectra of late-type galaxies show emission lines, so examining the spectra (e.g., using the PCA 
classification) is a simple way of removing such objects from the sample. Because the aperture of an SDSS 
spectroscopic fiber (3 arcsec) samples only the inner parts of nearby galaxies, and because the spectrum of 
the bulge of a nearby late-type galaxy can resemble that of an early-type galaxy, it is possible that some 
nearby late- type galaxies could be mistakenly included in the sample (e.g., Kochanek, Pahre & Falco 2000). 
Most of these will have been excluded by the first two cuts on the shape of the light profile. To check this, 
we visually inspected all galaxies with rdev > 8 arcsec. About ~ 50 of 225 (i.e., about 20%) looked like 
late-types, and so we removed them. Note that weak emission lines, such as H Q and/or O II, are still present 
in the early-type galaxy spectra in our sample. For example, the median H Q equivalent width for galaxies 
in our sample is —1.4 A, but ~ 5% of the galaxies have equivalent widths larger than 1 A and ~ 10% of 
these (i.e. 0.5% of the sample) have equivalent widths in the range 4 — 10 A. Similarly, although the median 
equivalent width of O II is 1.5 A, ~ 15% of the galaxies show O II equivalent widths larger than 4 A and 
~ 10% of these (i.e., 1.5% of the sample) have equivalent widths ranging from 8 to 15 A. 

The reason for our redshift cut is as follows. The SDSS main galaxy sample from which we select 
our early-type galaxies is apparent magnitude limited in r* , where the limits are defined using Pctrosian 
magnitudes m p . This means that the main galaxy sample is not magnitude limited in the other bands. 
Blanton et al. (2001) describe the cuts in m p one must apply in the other bands to obtain complete 
magnitude limited samples. They also note that m,j cv w 2.51og 10 0.8 + m p . Since we will almost always 
be working with m^ cv rather than m p , our cuts are slightly different from theirs; the cuts we use to define 
magnitude limited samples in the different bands are summarized in Table 1 (also see Figure 2). 

At z > 0.3, however, most early-type galaxies in the SDSS database were targeted using different 
selection criteria than were used for the main SDSS galaxy sample (see Strauss et al. 2002; they make-up 
the Luminous Red Galaxy sample described by Eisenstein et al. 2001). In the interest of keeping our sample 
as close to being magnitude limited as possible, we restricted our sample to z < 0.3. In addition, one might 
expect an increasing fraction of the early-type population at higher redshifts to have emission lines: if so, 
then our removal of emission line objects amounts to a small but redshift dependent selection effect. Since 
our sample is restricted to z < 0.3, this bias should be small. 

Before moving on, it is worth pointing out that there is a morphologically based selection cut which we 
could have made but didn't. Elliptical galaxies are expected to have axis ratios greater than about 0.6 (e.g., 
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Fig. 1. — Effective angular sizes rd ev , effective circular physical sizes R , velocity dispersions a, redshifts z, 
and (<?* — r*) colors as a function of axis ratio b/a for the galaxies in our r* sample. Bottom right panel 
shows that the typical axis ratio is b/a ~ 0.8. There is only a weak tendency for galaxies with small r^cv to 
be rounder, suggesting that the estimate of the shape is not compromised by seeing (typical seeing is about 
1.5 arcsec). Results in the other bands are similar. 
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Binney & Tremaine 1987). Since we have axis ratio measurements of all the objects in our sample, we could 
have included a cut on b/a. The bottom right panel of Figure 1 shows the distribution axis ratios b/a in our 
r* sample: about 20% of the objects in it have b/a < 0.6. (The spike at b/a = 1 is artificial — the pipeline 
reports large uncertainties in the fitted value of b/a for these objects.) 

Our combination of cuts on the shapes of the light profiles and spectral features mean that these objects 
are unlikely to be late-type galaxies. Indeed, a visual inspection of a random sample of the objects with 
axis ratios smaller than 0.6 shows that they look like SOs. The bottom left panel of Figure 1 shows that 
b/a does not correlate with color: in particular, the colors of the most flattened objects are not bluer than 
in the rest of the sample. Also, recall that all objects with r^cv > 8 arcsec were visually inspected and 
these, despite having b/a < 0.6 (top left panel), did not appear peculiar. In addition, b/a does not correlate 
with surface brightness or apparent magnitude. Galaxies with small angular sizes rdev are assigned large 
values of b/a only slightly more often than average (top left panel; the median b/a is 0.79, 0.78, 0.76, and 
0.7 for rd C v in the range 1-2, 2-3, 3-4 and greater than 4 arcsec). However, there is a weak trend for the 
objects at higher z to be rounder (middle right panel). This trend with z may be related to the magnitude 
limit of our sample rather than reflecting problems associated with the fits to the observed light profiles. 
(The two bands at z <~ 0.08 and z ~ 0.13 show that the effects of large scale clustering in our sample are 
significant.) Section 3.2.1 describes how we convert from the observed half-light radius r<j cv to an estimate of 
the physical half-light radius R a , and Section 3.3.1 discusses how we estimate velocity dispersions. Objects 
with smaller values of R tend to be more flattened (top right panel), and to have slightly smaller velocity 
dispersions (middle left); because of the magnitude limit, these objects drop out of our sample at higher 
redshifts. Requiring that b/a > 0.6 would remove such objects from our sample completely. 

In Paper III we study the Fundamental Plane populated by the galaxies in this sample. Excluding all 
objects with b/a < 0.6 has no effect on the shape of this Plane. So, in the interests of keeping our sample as 
close to being magnitude limited as possible, we chose not to make an additional selection cut on b/a. 

3.2. The photometric data of the sample 

3.2.1. Conversion to restframe luminosities and sizes 

Figure 2 shows the distribution of the m,d ev apparent magnitudes in our sample as a function of redshift. 
The dashed lines in each panel show the magnitude limits summarized in Table 1; the magnitude limited 
samples span the redshift range 0.01 < z < 0.3. 

To convert the apparent magnitude m to an absolute magnitude M we must assume a particular 
cosmology and account for the fact that at different redshifts an observed bandpass corresponds to different 
restframe bands (the K-correction). We write the Hubble constant today as lOO/i km s _1 Mpc _1 and use 
(Om, £^a, h) = (0.3, 0.7, 0.7). Most of our sample is at cz > 9000 km s _1 ; since line-of-sight peculiar velocities 
are not expected to exceed more than a few thousand km s -1 , we feel that it is reasonable to assume that all 
of a galaxy's redshift is due to the Hubble recession velocity. This means that we can compute the absolute 
magnitude in a given band by M = m — 51og 10 [Z?L(z; Om; ^a)] — 25 — K(z), where m is the apparent 
magnitude, is the luminosity distance in Mpc (from, e.g., Weinberg 1972; Hogg 1999), and K(z) is the 
K-correction for the band. 

Because we have five colors and a spectrum for each galaxy, we could, in principle, compute an empirical 
K-correction for each galaxy. This requires a good understanding of the accuracy of the SDSS photometry 
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Redshift 



Fig. 2. — Apparent magnitude (mdev) in the g* , r* , i* , and z* bands, versus redshift for galaxies which 
satisfied our selection criteria. In each band, magnitude limited samples are defined by including only those 
galaxies which fall between the dashed lines shown; these magnitude limits are summarized in Table 1. 
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Fig. 3. — Distribution of (seeing-corrected) effective angular sizes of galaxies in our sample. Typical seeing 
is about 1.5 arcsec. The distribution of effective radii in all the bands are very similar, although the radii 
are slightly larger in the bluer bands. 
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and spectroscopy, and should be possible when the survey is closer to completion. Rather than follow the 
procedure adopted by the 2dFGRS (Madgwick et al. 2002), or a procedure based on finding the closest 
template spectrum to each galaxy and using it to compute the K-correction (e.g., Lin ct al. 1999 for the 
CNOC2 survey), we use a single redshift dependent template spectrum to estimate the K-correction. In 
effect, although this allows galaxies at different redshift to be different, it ignores the fact that not all 
galaxies at the same redshift are alike. As a result, the absolute luminosities we compute are not as accurate 
as they could be, and this can introduce scatter in the various correlations we study below. Of course, using 
a realistic K-correction is important, because inaccuracies in K{z) can masquerade as evolutionary trends. 

For the reasons discussed more fully in Appendix A, our K-corrections are based on a combination of 
Bruzual & Chariot (2003) and Coleman, Wu & Wccdman (1980) prescriptions. Specifically, the K-corrections 
we apply were obtained by taking a Bruzual & Chariot model for a IO^Mq object which formed its stars 
with an IMF given by Kroupa (2001) in a single solar metallicity and abundance ratio burst 9 Gyr ago, 
computing the difference between the K-correction when evolution is allowed and ignored, and adding this 
difference to the K-corrections associated with Coleman, Wu & Weedman early type galaxy template. The 
results which follow are qualitatively similar for a number of other K-correction schemes (see Appendix A 
for details). 

In addition to correcting the observed apparent magnitudes to absolute magnitudes, we must also apply 
two corrections to convert the (seeing corrected) effective angular radii, rd ov , output by the SDSS pipeline 
to physical radii. First, we define the equivalent circular effective radius r Q = y/b/ardev- (Although the 
convention is to use r e to denote the effective radius, we feel that the notation r a is better, since it emphasizes 
that the radius is an effective circular, rather than elliptical aperture.) The reason we must make a second 
correction is shown in Figure 3. The different panels show the distribution of (seeing corrected) effective 
angular sizes rdev of the galaxies in the different bands. Notice that r^ev for most of the objects is larger than 
the typical seeing scale of 1.5 arcsec: we discuss this further in Section 3.2.2. Comparison of the mean sizes 
in the different panels (the text in the top right corner of each panel) shows that galaxies appear slightly 
larger in the bluer bands. 

Because our sample covers a reasonably large range in redshift, this trend means we must correct 
the effective sizes to a fixed restframe wavelength. Therefore, when converting from effective angular size 
r to effective physical size R we correct r (and the Petrosian radii rso and rgo) in each band by linearly 
interpolating from the observed bandpasses to the central rest wavelength of each filter. The typical correction 
is of the order of 4%, although it is sometimes as large as 10%. In this respect, this correction is analogous 
to the K-correction we would ideally have applied to the magnitude and surface brightness of each galaxy. 
Note that the dependence of size on wavelength is also important for analyses of color gradients and the 
color-magnitude relation presented in Paper IV. 

Our study will also require the effective surface brightness fi = — 2.51og 10 / o , where I is the mean 
surface brightness within the effective radius R (as opposed to the surface brightness at R ). In particular, 
we set /i Q = mdov + 2.5 log 10 (27rrg) — K(z) — 101og 10 (l + z). Note that this quantity is K-corrected, and 
also corrected for the cosmological (1 + z) 4 dimming. Our earlier remarks about the K-correction are also 
relevant here. 
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Fig. 4. — Distribution of ratio between the effective radius rd cv and tpsf, the FWHM of the PSF; the 
effective radius is typically about a factor of two larger than the seeing. 
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3.2.2. Reliability of photometric parameters 

There is some discussion in the literature about the effects of seeing on estimates of the effective radius, 
and hence the effective surface brightness. These effects are unimportant if the seeing scale, rpgp (the FWHM 
in arcsec), is an order of magnitude smaller than the effective radius (e.g., Saglia et al. 1993, 1997). Figure 4 
shows histograms of rd cv /fpsF in the g* , r*, i* and z* bands: the effective radius is typically only a factor 
of two larger than the seeing. Therefore, the accuracy of our estimates of the effective radii depend crucially 
on the correction for seeing being accurate. 

To illustrate that the de Vaucouleurs model fits output by the SDSS photometric pipeline are reasonable, 
Figure 5 shows the ratio of [HqcIa(z) / c] times the half-light radius to rpsF in the g* and r* bands (top and 
bottom panels) as a function of redshift, in a few volume limited catalogs. The catalog limits were chosen 
to be 0.5 mags wide, so the galaxies in any given catalog should all be quite similar to each other. (The 
factor H cIa(z)/c removes the expected scaling of angular diameter with redshift.) It is important to use 
volume limited samples to make the measurement, because size is known to correlate with luminosity (more 
luminous galaxies are larger on average), and in a magnitude limited sample such as ours, only the more 
luminous objects are present at higher redshifts. 

In the panels on the left, the half-light radius is estimated by rd cv , which has been corrected for seeing, 
whereas the panels on the right use r^, which is not seeing-corrected. Seeing tends to increase the half-light 
radius, and indeed, in the panel on the right, there is a clear trend for the half-light radii to increase with 
redshift. This trend is absent in the panels on the left, suggesting that the fitting procedure accounts quite 
well for the effects of seeing. The panels on the left even show a slight decrease with redshift we argue later 
(see discussion of Figure 12 in Section 4) that this may be a consequence of luminosity evolution. 

We have argued that our estimates of rd cv are probably not strongly compromised by seeing. In our 
dataset, rd ov is about a factor of two smaller than the radius r 10 % at which the fitted surface brightness 
is 10% of the sky level. Therefore, errors in the sky subtraction may also compromise our estimates of 
rdcv If the same object is moved to higher redshift but the sky level is kept fixed, then one expects 
7.67 (r w %/rdcv) 1 ^ 4 = const — 41n(l + z) because cosmological surface brightness dimming scales as (1 + z) 4 . 
Therefore, a plot of (riQ^Jr^cv) 1 ^ 4 versus ln(l + z), measured from volume limited catalog, should have a 
slope of —4/7.67. However, because each volume limited catalog spans a large range in redshift, we must also 
account for the surface brightness dimming associated with the fact that redshift dependent if-corrections 
must be made to estimate fixed restframe luminosities. Therefore, we expect the scaling with redshift to be 
slightly stronger: (r 10% /r dc v) 1/4 = const - 0.52 ln(l + z) - 0.052 ln(10) K(z). 

Figure 6 shows how (r 10 %/rdov) 1//4 + 0.12 K(z) depends on ln(l + z), in the same volume limited catalogs 
as before; solid lines show the mean trend in the faintest (top) to the most luminous (bottom) catalogs. 
Comparison with the dashed line, which has a slope of —0.52, shows that the scaling with redshift is close 
to that expected. (Note that it is important to use volume limited catalogs to make the measurement: the 
trend in any given volume limited subsample is shallower than the trend one would infer by first combining 
together all the subsamples and then making the measurement. To see why this happens, suppose that the 
mean size at fixed luminosity scales as R a oc L a . Then the magnitude limit induces a correlation which 
should scale approximately as rd ov <x Ro/z cx L a /z oc z 2a /z. If a > 1/2 then r^cv increases with z, so 
r io%/ r dov decreases with z simply because of the magnitude limit.) 

Some of the small differences in slope between the dashed and solid lines can be attributed to the effects 
of luminosity evolution. By fixing L, we select objects which are slightly smaller at higher z (c.f. Figure 12), 
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Fig. 5. — Effect of seeing on estimates of the half-light radius in the g* and r* bands (top and bottom panels). 
Each panel shows the ratio of [H dA{z)/c] times the angular half-light radius to tpsf as a function of redshift, 
in volume limited catalogs chosen to be 0.5 mags wide. The factor HqcIa(z) / c removes the expected scaling 
of angular diameter with redshift. In the panels on the left, the half-light radius shown is rd cv , which has 
been corrected for seeing, whereas the panels on the right use rso, which is not seeing-corrected. In the 
panels on the right, there is a clear trend for the half-light radii to increase with redshift, suggesting that 
seeing is affecting the estimates. This trend is absent in the panels on the left. (The panels on the left show 
a slight decrease with redshift — this may be a consequence of luminosity evolution.) 
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Fig. 6. — The ratio of the radius r 10 % s ky at which the fitted surface brightness is 10% of the sky level, and 
the effective radius rdcv, as a function of redshift, for a few volume limited subsamples. In all four bands, 
this ratio decreases with increasing redshift, although the trend is slightly weaker in g* than in the other 
bands. Dashed line shows the expected scaling if there is no luminosity or size evolution. The trend for 
the combination of all the samples is substantially steeper than that within any one volume limited sample, 
suggesting that the typical size increases strongly with luminosity. 
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hence have slightly smaller values of rdev and so larger values of rio%/rdcv than expected; so they appear to 
decline less rapidly than const— 0.52 z. Luminosity evolution is expected to be stronger in the bluer bands, 
and the g* band relation does appear to be slightly shallower than the other bands. Thus, as was the case 
for rdov/^PSF, we see redshift dependent trends in r 10 %/rd C v which we believe are physically reasonable. The 
sky-subtraction procedure does not appear to have introduced any obvious artifacts into our estimates of 

Finally, we turn to an inspection of how well the de Vaucouleurs model actually fits the observed surface 
brightness profiles. Figure 7 shows angular averaged surface brightness profiles in the g* , r* , i* and z* bands 
(bottom to top in each of the upper panels), for a random selection of galaxies in our sample, ordered 
by redshift and size. The text in each panel shows the plate number, MJD, and fiber ID number of the 
observation, the redshift z, the fitted apparent r* magnitude m r * , the fitted axis ratio b/a, and the effective 
physical radius R , computed from the fitted angular radius rdev, the fitted axis ratio and the redshift, as 
described previously. The solid lines show de Vaucouleurs profiles which, when convolved with the pipeline's 
model for the seeing, provide the best fit to the symbols (because the fits are actually not to these azimuthally 
averaged one-dimensional profiles, but to the two dimensional profile, they serve mainly as guides to the eye.) 
The horizontal lines show the surface brightness at 1% of the sky in the different bands. Long, medium and 
short tick marks on these lines show rpsF, rd ov and r w %, the radius at which the fitted surface brightness 
equals 10% of the sky. 

The bottom panels show the residuals (observed minus fit) in the different bands (crosses, stars, dia- 
monds and triangles show the g* , r* , i* and z* band results, respectively). The residuals are systematically 
low on small scales, but comparison with the large tick marks, which indicate the size of the seeing disk, 
show that this discrepancy is because we are comparing fits which have been corrected for seeing with the 
observed surface brightnesses which are not seeing-corrected. On larger scales (typically larger than twice the 
half-light radii), the residuals are also large, indicating that there may be problems with the sky subtraction. 
In linear flux units, the residuals tend to a constant value which is typically smaller than 0.3% of the sky 
level (recall that the horizontal lines in the panels show the surface brightness at 1% of sky). Although 
this is not an entirely fair comparison, since the fits are not to performed on these azimuthally averaged 
one-dimensional profiles, but to the two dimensional profile, this test suggests that the fit is reliable, but 
that the zero-point level of the sky may be systematically biased by ~ 0.3%. Note that the residuals are 
larger in the g* and z* bands; but this is consistent with the larger errors we quote on the photometry in 
these bands. 

Figure 8 shows the r* images of the same objects whose profiles are shown in Figure 7. The images are 
each 11.88 x 11.88 arcsec 2 (i.e., 30 pixels on a side). One might expect some of the worst fits are due to the 
presence of a disk component which our de Vaucouleur model does not account for. However, comparison 
of the images and the fits shows no obvious correlation between goodness-of-fit and, e.g., b/a. 

In summary, we have performed a number of tests of the SDSS photometric pipeline reductions and have 
found little evidence that they are significantly biased. Therefore, in what follows, the apparent magnitudes 
we use are those output by the SDSS photometric pipeline. To convert to luminosities, we K-correct (and 
use the inverse square law). The angular sizes are output from the SDSS photometric pipeline. We correct 
these to a fixed restframe size, and we correct for a/b. Both these corrections are standard. Thereafter, 
surface brightness is defined from the restframe luminosity and size with no additional corrections. The 
uncertainties on the photometric parameters we use in this paper, and the magnitude limits in each band, 
are summarized in Table 1. There is some subtlety in obtaining them from the format available in the SDSS 
database; this is discussed in Appendix D. 
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Fig. 7. — Examples of surface brightness profiles selected randomly from galaxies in our sample which are 
in the redshift range z < 0.075. (Upper panels) Text in each panel shows plate number, MJD, and fiber 
ID of the object, the redshift z, the apparent r* magnitude m r », the fitted axis ratio b/a, and the physical 
half-light radius R . Panels are ordered by physical radius R 0l which increases from top to bottom, and 
from left to right. Symbols show the surface brightness profiles in the g* , r* , i* and z* bands (bottom to 
top in each of the upper panels), and solid lines show the de Vaucouleur fits which, when convolved with 
the seeing, provide the best fit to the symbols. Horizontal lines show 1% of the surface brightness of the 
sky in the different bands, and large, medium and small tickmarks show rpgF, r<iev an d fio%, the radius at 
which the fitted surface brightness equals 10% of the sky. (Bottom panels) The residuals (observed minus 




Fig. 8. — r*-band images for the z < 0.075 objects shown in Figure 7. The labels on the x- and y-axes show 
the coordinates (RA and DEC in J2000) of the galaxy centers. 
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Fig. 7. - Continued. Redshift range 0.075 < z < 0.1. 
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Fig. 8. - Continued. Redshift range 0.075 < z < 0.1. 



-22 - 



0.1 < z < 0.12 




Fig. 7. - Continued. Redshift range 0.1 < z < 0.12. 



- 23 - 




-24- 



0.12 < z < 0.14 




'IIS 51900 -108 
z = 0.133 
m r .= 1 7.4 
b/a = 0.90 
R„= 4.5 




440 51885 558 
z = 0.127 
m r .= 1 6.8 
b/a = 0.98 
R„= 5.4 




■119 51868 15 
z = 0.132 
m r .= 1 6.8 
b/a = 0.95 
R„= 4.8 




449 51900 451 

z = 0.123 
m r .= 16.8 
b/a = 0.90 
R = 5.5 




291 51660 62 
Z = 0.126 
m r .= 17.1 
b/a = 0.95 
R„= 5.4 





I i\ * 




442 51882 226 
z = 0.12b 
m r .= 1 6.4 
b/a = 0.77 
R„= 7.0 




305 51613 402 
z = 0.134 
m r .= 16.3 
b/a = 0.95 
R„= 8.1 



2.0 1.0 1.5 

r Foresee! 7/4 



Fig. 7. - Continued. Redshift range 0.12 < z < 0.14. 
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Fig. 7. - Continued. Redshift range 0.14 < z < 0.18. 
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Fig. 7. - Continued. Redshift range z > 0.18. 
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Continued. Redshift range z > 0.18. 
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Table 1: Photometric parameters and median errors of the objects in our sample. 
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3.3. The spectroscopic data of the sample 

3.3.1. Estimating the velocity dispersions 

The SDSS pipeline does not provide an estimate of the line-of-sight velocity dispersion, a, within a 
galaxy, so we must compute it separately for the early-type galaxy sample. The observed velocity dispersion 
a is the result of the superposition of many individual stellar spectra, each of which has been Doppler 
shifted because of the star's motion within the galaxy. Therefore, it can be determined by analyzing the 
integrated spectrum of the whole galaxy. A number of objective and accurate methods for making velocity 
dispersion measurements have been developed (Sargent et al. 1977; Tonry & Davis 1979; Franx, Illingworth 
& Heckman 1989; Bender 1990; Rix & White 1992). Each of these methods has its own strengths, weaknesses, 
and biases. Appendix B describes how we combined these different techniques to estimate a for the galaxies 
in our sample. 

The velocity dispersion estimates we use in what follows are obtained by fitting the rcstframe wavelength 
range 4000 - 7000 A, and then averaging the estimates provided by the Fourier-fitting and direct-fitting 
methods to define what we call cr cs t. (We do not use the cross-correlation estimate because of its behavior 
at low S/N as discussed in Appendix B. Note that the S/N of the SDSS spectra depend on wavelength.) 
The error on cr cs t is determined by adding in quadrature the errors on the two estimates (i.e., the Fourier- 
fitting and direct-fitting) which we averaged. The resulting error is between <51ogcr <~ 0.02 dex and 0.06 dex, 
depending on the signal-to-noise of the spectra, with a median value of ~ 0.03 dex. A few galaxies in our 
sample have been observed more than once. The scatter between different measurements is ~ 0.04 dex, 
consistent with the amplitude of the errors on the measurements (see Figure 19). Based on the typical S/N 
of the SDSS spectra and the instrumental resolution, we chose 70 km s _1 as a lower limit on the velocity 
dispersions we use in this paper. 

Following J0rgensen et al. (1995) and Wegner et al. (1999), we correct a cs t to a standard relative 
circular aperture defined to be one-eighth of the effective radius: 

/ x 0.04 

^cor / ^fiber \ /-. \ 

^ ~ \7j8 J ' 1 ' 

where rgbcr = 1-5 arcsec and r Q is the effective radius of the galaxy measured in arcseconds. In principle, 
we should also account for the effects of seeing on cr cs t, just as we do for r a . However, because the aperture 
correction depends so weakly on r a (as the 0.04 power), this is not likely to be a significant effect. In any 
case, most galaxies in our sample have r Q > 1.5 arcsecs (see Figure 3). 
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Note that this correction assumes that the velocity dispersion profiles of early-type galaxies having 
different r a are similar. At the present time, we do not have measurements of the profiles of any of the 
galaxies in our sample, so we cannot test this assumption. In Paper II of this series we will argue that the 
galaxies in our sample evolve very little; this means that if we select galaxies of the same luminosity and 
effective radius, then a plot of velocity dispersion versus redshift of these objects should allow us to determine 
if the aperture correction above is accurate. The results of this exercise are presented in Appendix C. 

3.3.2. Examples of early-type galaxy spectra 

To illustrate the typical quality of the spectra in our sample, Figure 9 shows spectra of the galaxies 
which lie along the diagonals in Figure 7 (panels from top to bottom here are for objects which are top 
left to bottom right in figure 7). All the spectra in our sample have PCA classification numbers a < —0.1, 
typical of early- type galaxy spectra (Connolly & Szalay 1999), although a few of the spectra do show weak 
emission lines. 

The typical signal-to-noise ratio of the spectra in our sample is 16, although S/N does vary with 
wavelength. Therefore, individual spectra cannot be used to make reliable estimates of, e.g., the Lick 
indices. However, by co-adding spectra of similar objects (i.e., galaxies with similar luminosities, sizes, 
velocity dispersions, and redshifts) it is possible to construct composite spectra which have considerably 
better S/N. This is the subject of Paper IV. 

3.4. An estimate of the local density 

In subsequent papers, we will be interested in how the properties of the galaxies in our sample depend 
on their local environment. To do so, we must come up with a working definition of environment. The 
set of galaxies in the SDSS photometric database is much larger than those for which the survey actually 
measures redshifts. Some of these galaxies may well lie close to galaxies in our sample, in which case they 
will contribute to the local density. We would like to find some way of accounting for such objects when we 
estimate the local density. 

For a subset of the galaxies in our sample, the colors expected of a passively evolving early-type were used 
to select a region in g*-r* versus r*-i* color space at the redshift of the galaxy of interest. All galaxies within 
0.1 magnitudes in color of this point were included if they were: a) within XhT 1 Mpc of the main galaxy, and 
b) brighter than —20.25 in Mi-. (The box in color space is sufficiently large that the difference between this 
techinique, and using the observed colors themselves to define the selection box is not important.) These 
two cuts are made assuming every galaxy in the color-color range is at the redshift of the galaxy of interest. 
The end result of this is that each galaxy in the subsample is assigned a number of neighbors. Note that, 
because of the selection on color, our estimate of the local density is actually an estimate of the number 
of neighbors which have the same colors as early-type galaxies. In the papers which follow, we will often 
present results for different redshift bins. When we do, it is important to bear in mind that this procedure 
for assigning neighbours is least secure in the lowest redshift bin (typically z < 0.08). 

Figure 10 shows how the luminosities, surface brightnesses, sizes, velocity dispersions and (a combination 
of the) axis ratios depend on environment. The different symbols for each bin in density show averages over 
galaxies in different redshift bins: circles, diamonds, triangles, and squares are for galaxies with redshifts in 
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Fig. 9. — Spectra of galaxies which lie along the diagonals of Figures 7 and 8. Text in each panel gives the 
plate number, MJD and fiber ID, the redshift z, the signal-to-noise ratio, and the PCA classification number. 
Some common absoption features are also indicated. The line along the bottom of each panel shows the 
error in the flux density in each pixel. Here, galaxies span the redshift range z < 0.075. 




Fig. 9. - Continued. Redshift range 0.075 < z < 0.1. 




Fig. 9. - Continued. Redshift range O.K z < 0.12. 




Fig. 9. - Continued. Redshift range 0.12 < z < 0.14. 




Fig. 9. - Continued. Redshift range 0.14 < z < 0.18. 




Fig. 9. - Continued. Redshift range z > 0.18. 
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the range 0.075 < z < 0.1, 0.1 < z < 0.12, 0.12 < z < 0.14, and 0.14 < z < 0.18. Error bars show the error 
in the mean value for each bin. Symbols for the higher redshift bins have been offset slightly to the right. 

For any given set of symbols, the bottom right panel shows that the mean redshift in each bin in density 
is not very different from the mean redshift averaged over all bins. This suggests that our procedure for 
estimating the local densities is not biased. The other panels show corresponding plots for the other observed 
parameters. When the number of near neighbours is small, the luminosities, sizes and velocity dispersions 
all increase slightly as the local density increases, whereas the surface brightnesses decrease slightly. All 
these trends are very weak. The bottom right panel shows y/a/b — 1— < \J a/b — 1 >, where b/a is the 
axis ratio. [Some authors (e.g. dc Zccuw & Franx 1991) suggest that this is a good approximation of the 
ratio of rotational to random motions within the galaxy, (v/a)i so , for oblate, rotationally flattened galaxies. 
Bright ellipticals have been shown to have a mean (v/a)/(v/<j)i so ~ O.4.] There are no obvious trends with 
environment. It is difficult to say with certainty that the trends with environment in the top four panels of 
Figure 10 are significant. A more efficient way of seeing if the properties of galaxies depend on environment 
is to show the residuals from the Fundamental Plane. This is done in Paper III. 

Figure 11 shows a redshift-space pie-diagram distribution of our sample. Most of the sample is at low 
declination (\S\ < 2°); in addition, there are three wedges from three different disconnected regions on the 
sky. Red and blue symbols denote galaxies which were classified as being in dense and underdense regions 
(as described in Section 3.4 below), whereas black symbols show galaxies in groups of intermediate richness, 
or for which the local density was not determined. 

4. The early-type galaxy catalog 

The observed parameters of each galaxy in our early-type galaxy sample are given in Table 4. The 
different columns show the ra, dec, redshift z and its error estimate, the signal-to-noise ratio S/N of the 
spectrum, the apparent half light radius r<j ov and error estimate, the fitted apparent magnitude rtidev and 
error estimate, the model magnitude m m and error estimate, and the axis ratio b/a and error estimate. The 
photometric parameters in this table are those measured in the r* band. The complete version of this table, 
and similar tables of the catalog parameters in the g* , i* and z* bands are available in the electronic edition 
of the Journal. Because we will be interested in the distribution of colors for a galaxy sample defined in, e.g., 
the r* band, the tables present the parameters of the same galaxies in all the bands. To select the magnitude 
limited samples we use in subsequent papers refer to Table 1. The errors (5mdcv- P hoto and <5m mo d- p hoto are 
the error estimates output by the SDSS photometric pipeline. As discussed in Appendix D, they are not 
appropriate estimates of the errors in TOdcv and m mo dei- O ur procedure for making more appropriate error 
estimates is described in Appendix D. 

Table 4 shows restframe quantities we computed from the quantities in Table 4. The different columns 
show the ra, dec, aperture corrected velocity dispersion log 10 a and its error estimate, the physical size 
log 10 R and error estimate, the surface brightness \x and error estimate, the absolute magnitude Md ev and 
error estimate, the absolute model magnitude M mo d and error estimate, the K-correction, and the number 
of neighbours which we use as an estimate of the local density. As for Table 4, a complete version of this 
table is in the electronic edition of the Journal, as are similar tables showing parameters in the g* , i* and z* 
bands. (The errors in Md cv and M mo d c i listed in this table were derived from the estimates <5m p hoto of the 
previous table. As discussed in Appendix D, they are also the appropriate estimates for the error on md cv 
and m mod .) 



- 39 - 





0.4 








0.2 


A 




N 







-0.2 


V 


-0.4 


A 




N 


0.04 



b 

V 



b 

o 



too 

O " 
J A 



V 



- zM, 



35 





0.04 
0.1 

0.1 



"i r 




4 



"3, <> 



J L 



■ > 



"i r 



i i 




1 



[] - 



J L 



0.4 
0.2 



0.2 
0.4 



0.1 



-0.1 
0.01 


-0.01 



A 

N 

V 

I 



tr -1 
o 

o 



A 

N 

V 

N 

I 

A 

N 

V 



20 40 
Ndensity 



20 40 
Ndensity 



Fig. 10. — Luminosities, surface brightnesses, velocity dispersions, sizes, axis ratios, and mean redshifts, as 
a function of nearby early- type neighbours. The different symbols for each bin in density show averages over 
galaxies in different redshift bins: circles, diamonds, triangles, and squares are for galaxies with redshifts in 
the range 0.075 < z < 0.1, 0.1 < z < 0.12, 0.12 < z < 0.14, and 0.14 < z < 0.18. Although the velocity 
dispersions appear to increase with increasing local density, the increase is small. 
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Fig. 11. — Pie-diagram distribution of our sample. Most of the sample is at low declination (|<5| < 2°), but 
three wedges are at higher declinations (as indicated). Right ascension increases clockwise, with the zero 
at twelve o'clock. Galaxies with many (> 15) and a few (< 2) near neighbours are shown with red and 
blue dots, whereas those in the intermediate regime, or for which the local number of neighbours was not 
determined, are shown with black dots. 



Tabic 2. The SDSS early-type galaxies catalog: Observed parameters 



RA (J2000) 
(deg) 


DEC (J2000) 

(deg) 




z 




Sz 


S/N 


1 dev 


u 1 dev 


777 ,1^,.. 

l^magj 


f) T77 ,i i j. „ 

u it ^dev — photo 

1 m o cy 1 

l^magj 


77 7 ™ . . , j 

'"'mod 

I YlTi Q <T 1 

l^magj 


i) 777 „ j i j. „ 

u " t moa-pnoto 


b/a 


5b /a 


113.810738 


36.307522 


0. 


.25100 


0. 


,00009 


12 


3.45 


0.15 


17.234 


0.009 


17.234 


0.009 


0.82 


0.03 


115.321098 


37.511093 


0. 


.08513 


0. 


.00008 


19 


2.55 


0.04 


16.400 


0.009 


16.400 


0.009 


0.99 


0.01 


114.201271 


36.971313 


0. 


.07832 


0. 


.00008 


15 


2.16 


0.05 


17.202 


0.009 


17.202 


0.009 


0.71 


0.02 


114.306282 


37.111927 


0. 


.13049 





,00008 


20 


2.48 


0.05 


16.828 


0.008 


16.828 


0.008 


0.68 


0.01 


113.294510 


36.370544 


0. 


.21756 


0. 


.00007 


13 


2.18 


0.08 


17.451 


0.008 


17.451 


0.008 


0.80 


0.03 


113.376007 


36.478168 


0. 


13365 


0. 


.00008 


14 


2.45 


0.08 


17.304 


0.008 


17.304 


0.008 


0.63 


0.02 


113.335358 


36.648064 


0. 


.14016 


0. 


,00007 


21 


2.01 


0.04 


16.871 


0.006 


16.871 


0.006 


0.68 


0.01 


113.902443 


37.023502 


0. 


.20225 


0. 


,00008 


12 


4.17 


0.17 


16.998 


0.011 


16.998 


0.011 


0.96 


0.03 


114.022346 


37.099522 


0. 


.17845 


0. 


,00008 


14 


2.07 


0.08 


17.453 


0.011 


17.453 


0.011 


0.76 


0.03 


113.925407 


37.113167 


0. 


.18365 


0. 


,00009 


15 


2.71 


0.08 


17.053 


0.010 


17.053 


0.010 


0.80 


0.02 


114.060677 


37.211426 


0. 


09610 


0. 


,00008 


13 


2.11 


0.07 


17.279 


0.010 


17.279 


0.010 


0.85 


0.03 


113.827904 


37.380821 


0. 


.18387 


0. 


,00008 


12 


4.82 


0.18 


16.898 


0.009 


16.898 


0.009 


0.97 


0.03 


113.945900 


37.708260 


0. 


.07476 


0. 


,00006 


19 


2.98 


0.05 


16.882 


0.009 


16.882 


0.009 


0.34 


0.01 



Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. This table Jists 
the r* band photometric parameters. Similar tables showing parameters in the g* , i* and z* bands are also available in the electronic edition. ^ 



Table 3. The SDSS early-type galaxies catalog: Physical parameters 



RA 


DEC 


V 


sv 


R 


SR 




5ji 


M dcv 


6M dcv 


M mod 


SM mod 


Kcorr 


Nden 


(deg) 


fdeg) 


(dex) 


(dex) 


(dex) 


(dex) 


(mag) 

V 0/ 


(mag) 


(mag) 


(mag) 


(mas) 


(mas) 

V 0/ 


(mag) 




113.810738 


36.307522 


2.501 


0.034 


0.979 


0.020 


19.873 


0.060 


-23.587 


0.033 


-23.587 


0.033 


0.309 


4 


115.321098 


37.511093 


2.133 


0.029 


0.600 


0.007 


19.918 


0.020 


-21.648 


0.018 


-21.648 


0.018 


0.104 


4 


114.201271 


36.971313 


2.055 


0.045 


0.414 


0.012 


19.993 


0.034 


-20.643 


0.021 


-20.643 


0.021 


0.093 


3 


114.306282 


37.111927 


2.363 


0.022 


0.664 


0.010 


19.620 


0.030 


-22.268 


0.019 


-22.268 


0.019 


0.159 


2 


113.294510 


36.370544 


2.307 


0.040 


0.871 


0.017 


19.944 


0.052 


-22.978 


0.029 


-22.978 


0.029 


0.268 


4 


113.376007 


36.478168 


2.319 


0.026 


0.619 


0.016 


19.810 


0.046 


-21.851 


0.026 


-21.851 


0.026 


0.163 


2 


113.335358 


36.648064 


2.266 


0.024 


0.541 


0.010 


18.866 


0.028 


-22.407 


0.017 


-22.407 


0.017 


0.173 


5 


113.902443 


37.023502 


2.352 


0.045 


1.110 


0.019 


20.886 


0.059 


-23.229 


0.034 


-23.229 


0.034 


0.244 


999 


114.022346 


37.099522 


2.388 


0.038 


0.749 


0.018 


19.872 


0.054 


-22.440 


0.031 


-22.440 


0.031 


0.212 


7 


113.925407 


37.113167 


2.433 


0.029 


0.856 


0.014 


19.934 


0.041 


-22.914 


0.025 


-22.914 


0.025 


0.217 


6 


114.060677 


37.211426 


2.173 


0.034 


0.542 


0.015 


20.219 


0.044 


-21.058 


0.027 


-21.058 


0.027 


0.114 


1 


113.827904 


37.380821 


2.347 


0.040 


1.094 


0.018 


20.965 


0.055 


-23.071 


0.030 


-23.072 


0.030 


0.217 


4 


113.945900 


37.708260 


2.090 


0.028 


0.394 


0.009 


19.687 


0.024 


-20.851 


0.019 


-20.851 


0.019 


0.088 


3 



Note. — RA and DEC are in J2000. We have set V — logger with a in km s R — log^o R a with R Q in kpc h^^, Kcorr is the 
K-correction and Nden is the number of near neighbors which we use to estimate the local density. Objects for which the number of 
neighbors was not determined have been assigned Ndcn=999, and we have set Nden=100 even if the number of neighbors is larger than 
100. The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. This table 
lists the r* band photometric parameters. Similar tables showing parameters in the g* , i* and z* bands are also available in the electronic 
edition. 
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Fig. 12. — Surface brightnesses, velocity dispersions, sizes, masses and densities of galaxies as a function of 
redsliift, for a few bins in r* luminosity. Top left panel shows volume limited catalogs which do not overlap 
in r* luminosity, dots in the other panels show the galaxies in the volume-limited subsamples defined by the 
top-left panel, and solid lines show the mean trend with redshift in each subsample. Results in g* , i* and 
z* are similar. 
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Some properties of our r* sample as a function of redshift z are shown in Figure 12. The panels show 
the K-corrected absolute magnitude M, the K-corrected effective surface brightness fj, , also corrected for 
cosmological surface brightness dimming, the effective circular radius R in h^Q kpc, corrected to a standard 
restframe wavelength, the aperture corrected velocity dispersion a in km s _1 , and two quantities which are 
related to an effective mass and density, all plotted as a function of redshift. 

The bold lines in the top left panel show the effect of the apparent magnitude cuts. There is, in addition, 
a cut at small velocity dispersion (<~ 70 km s _1 ) which, for our purposes here, is mostly irrelevant. The 
apparent magnitude cuts imply complex z-dependent cuts on the other parameters we observe. In the papers 
which follow, we attempt to account correctly for the selection effects that the magnitude cuts imply. 

For small intervals in luminosity, our sample is complete over a reasonably large range in redshifts. 
To illustrate, the thin boxes in the top-left panel show bins in absolute magnitude of width 0.5 mags over 
which the sample is complete. The solid lines in the other panels show how the median surface-brightnesses, 
sizes, velocity dispersions, masses and densities of galaxies in each r* luminosity bin change as a function 
of redshift. Although all these quantities depend on luminosity, the figure shows that, at fixed luminosity, 
there is some evidence for evolution: at fixed luminosity the average surface brightness is brightening. The 
size at fixed luminosity decreases at a rate which is about five times smaller than the rate of change of \x a . 
This suggests that it is the luminosities which are changing and not the sizes. (To see why, suppose that 
the average size at fixed absolute magnitude is (log 10 [i? o /i?»(z)]) = s[M — M*(z)], where i?* and M* arc 
characteristic values, and s is the slope of this mean relation. If the characteristic luminosity increases with 
z, but the characteristic size remains constant, R*(z) = R*(0), and the slope of the relation also does not 
change, then the mean size at fixed M decreases with z. The surface brightness is \i Q oc M + 5 log 10 R Q , so the 
mean \i a at fixed M changes five times faster than the mean R Q , at fixed luminosity. This argument is also 
relevant to Figures 5 and 6, since those plots use rd cv which differs from R Q by two small correction factors 
[c.f. Section 3.2.1].) In Paper II we argue that this trend is qualitatively consistent with that expected of a 
passively evolving population. 



5. Summary 

This is the first of four papers which study the properties of early- type galaxies at relatively low redshifts 
(z < 0.3). This paper describes how we selected the sample from the SDSS database using objective 
criteria and present the measured photometric and spectroscopic parameters for ~ 9000 early-type galaxies 
(Table 4). The database contains the redshift of these galaxies as well as reliable measurements of a number 
of photometric properties (luminosities, sizes, surface brightnesses, colors) in the g* , r* , i* and z* bands. 
(Data in the u* band is also available, but, because it is noisier, we did not use it here.) 

The galaxies in our sample span the redshift range 0.01 < z < 0.3. Therefore, differences between 
observed and restframe wavelengths are not negligible. We discussed a number of methods for estimating 
the appropriate K-correction which must be applied to obtain luminosities which sample the same restframe 
wavelength range (Appendix A). Since the half-light radii also depend on wavelength (Figure 3), we also 
corrected these to sizes at fixed restframe wavelengths. Finally, we estimated the internal velocity dispersions 
of these galaxies (Appendix B), which, following common practice, we corrected to an aperture of r e /8 (c.f. 
equation 1). The resulting catalog listing the restframe photometric and spectroscopic parameters of ~ 9000 
early-type galaxies (Table 4) is useful for measuring correlations between these various observables, such as 
the Faber-Jackson L — <r relation, the Kormendy R a — I a relation, the L — R D relation, the Fundamental 
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Plane R — I Q — a relation, and the color-magnitude relation. 

The sample is essentially magnitude limited (Table 1, Figure 2), and the galaxies in it populate a range 
of environments (Figure 11). Therefore, the sample can be used for studying how the correlations discussed 
above depend on environment and redshift. For instance, Figure 12 shows evidence for weak evolution in 
the luminosities, and Figure 10 shows evidence for some weak environmental dependences. These will be 
quantified in Papers II and III. 

Neither the K-corrections nor the aperture corrections we apply to the velocity dispersions are ideal. By 
the time the Sloan Digital Sky Survey is complete, the uncertainties in the K-corrections, which prevent us 
at the present time from making precise quantitative statements about the evolution of the luminosities and 
colors, will be better understood. In addition, the size of the sample will have increased by more than an 
order of magnitude. This larger sample will allow us to use the data directly to estimate the correct aperture 
correction which should be applied to the velocity dispersions. This is discussed more fully in Appendix C. 
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A. K-corrections 

When converting the observed apparent magnitude to the rest-frame absolute magnitude of an object, we 
must account for the fact that the SDSS filters measure the light from a fixed spectral range in the observers 
rest-frame; therefore, they measure different parts of the rest-frame spectrum of galaxies at different redshifts. 
Correcting for this is known as the K-correction. One way to make this correction is to assume that all the 
galaxies at a given redshift are similar, and to use an empirically determined template spectrum, measured 
from a few accurately measured spectra, to estimate the K-correction. Using a mean color to estimate the 
K-correction is not ideal. When the survey is closer to completion it should become possible to make this 
correction on an object-by-object basis. 

Empirically determined template spectra for early-type galaxies at low redshifts exist (e.g. Coleman, Wu 
& Weedeman 1980; Fukugita, Shimasaku & Ichikawa 1995). (We used N. Benitez's Bayesian Photometric 
Redshift package (Benitez 2000) to derive K-corrections for the Coleman, Wu & Weedman early-type galaxy 
template in the SDSS passbands.) If we were certain that evolution effects were not important, then these 
empirically determined K-corrections would allow us to work out the K-corrections we should apply to the 
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Fig. 13. — Difference between K-corrections based on two models (Coleman, Wu & Weedman 1980 and 
Bruzual & Chariot 2003) of the SDSS colors of early-type galaxies. Filled circles, crosses, stars and diamonds 
are for the g*, r*, i* and z* bands. 

high redshift population. However, if the stars in early-type galaxies formed at approximately the same 
time, and if the mass in the galaxies has remained constant, so the evolution is entirely due to the passive 
aging of the stellar population, then, the mass to light ratio of early-type galaxies is expected to vary 
approximately as M/L oc (t — tform) -0 ' 6 ( e -g-i Tinsley & Gunn 1976), with the precise scaling being different 
in different bands. Because our sample spans a relatively large range in redshift, we may be sensitive to 
the effects of this passive evolution. In addition, because the sample is large, we may be able to measure, 
and hence be sensitive to, even a relatively small amount of evolution. For this reason, it would be nice to 
have a prescription for making K-corrections which accounts for evolution. Absent empirically determined 
templates for this evolution, we must use stellar population systhesis models to estimate this evolution, and 
so determine K-corrections for different bands. 

As a first example, we chose a Bruzual & Chariot (2003) model for a 10 n Mo object which formed its 
stars with the universal IMF given by Kroupa (2001) in a single solar metallicity and abundance ratio burst 
9 Gyr ago. We then recorded how its colors, as observed through the SDSS filters, changed as it was moved 
through redshift without altering its age. This provides what we will call the no-evolution K-correction. 
Figure 13 shows a comparison of this with the empirical Coleman, Wu & Weedman (1980) nonevolving K- 
corrections. The two estimates are in good agreement in g* and i* out to z ~ 0.3. They differ substantially 
at higher redshifts, but this is not a concern because none of the galaxies in our sample are so distant. In r* 
and z* the two estimates agree only at z < 0.15. Therefore, quantitative estimates of evolution in luminosity 
and/or color will depend on which K-correction we use. 

Figure 14 compares both sets of nonevolving templates with the observed colors of the galaxies in our 
sample. The upper set of curves in each panel show the colors associated with the nonevolving Coleman, 
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Fig. 14. — Apparent colors of the galaxies in our sample. In each panel, dotted and solid lines show the non- 
evolving and evolving Coleman et al. templates, whereas dashed and dot-dashed show Bruzual & Chariot 
models. The upper set of curves in each panel show what one expects to see if the intrinsic colors of galaxies 
at higher redshifts are the same as they are nearby, whereas the lower sets of curves show the predictions if 
the higher redshift population is slightly younger, and so bluer. The magnitude limit of the sample makes 
it appear as though the no-evolution curves describe our data well. In the main text, we use the lower solid 
curve to make K-corrections to the observed magnitudes. 
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Wu & Weedman (1980) template (dotted) and the Bruzual & Chariot (2003) no-evolution model (dashed). 
The figure shows that both predictions for g* — r* are similar, but that they are different for r* — i* and 
r* — z*, with the differences increasing with redshift. [In both cases, we have shifted the predicted g* — r* 
blueward by 0.08 mag at all z. Such an offset appears to be required for the SDSS photometric calibrations 
in Stoughton ct al. (2002) which we use here (also see Eiscnstein ct al. 2001 and Strauss et al. 2002), 
although the reason for it is not understood.] 

We argue in Papers II, III and IV that the luminosities, colors and chemical abundance ratios of these 
objects show evidence for passive evolution: the higher redshift population appears to be slightly younger. 
Therefore, our next step is to include the effects of evolution. Because the predicted observed colors at 
redshift zero are in good agreement with our data, we took the same Bruzual & Chariot (2003) model, 
but this time we recorded how its rest-frame colors evolve with redshift, and then computed what these 
evolved (i.e., z-dependent) colors look like when observed in the SDSS filters. These evolving colors are 
shown as the dot-dashed lines in Figure 14. In an attempt to include evolution in the CWW templates, 
we set K^$ w (z) = K^ v w (z) + K<gg l (z) - K%°^(z). The lower solid lines in each panel of Figure 14 
show the observed g* — r* , r* — i* and r* — z* colors associated with these evolving models (and again, 
the predicted g* — r* curves have been shifted blueward (downward) by 0.08 mag at all z). Comparing the 
evolving Bruzual-Charlot and the Coleman et al. colors with the upper set of no-evolution curves shows the 
evolution towards the blue at high redshift. Although the data appear to be very well fit by the no-evolution 
curves, this agreement is slightly misleading. More luminous galaxies tend to be redder. As a consequence, a 
magnitude limited catalog contains only the redder objects of the higher redshift population. A curve which 
describes the colors of the population as a whole will therefore appear to be biased blue. 

Figure 15 shows this explicitly. The two panels were constructed by making mock catalogs of a passively 
evolving population (i.e., the higher redshift population is brighter and bluer) in which our estimates of 
the correlation between colors, luminosities, velocity dispersions and color and luminosity evolution were 
included (see Appendix A in Paper II) . The top panel shows the distribution of observed colors if there were 
no magnitude limit, and the bottom panel shows the observed colors of a magnitude limited sample. The 
solid curve, the same in both panels, is the predicted trend of color with redshift which we use to make our 
K-corrections; i.e., K^^ w (z). Notice that although it describes the complete simulations well, it is bluer 
than the higher redshift galaxies in the magnitude limited sample. Comparison with the previous figure 
shows that the difference here is similar to that seen in the real data, suggesting that our K-corrections and 
evolution estimates of the mean of the population are self-consistent. 

Of course, if we do not observe the mean of the high redshift population, but only the redder fraction, 
then we must decide whether it is realistic to use a K-correction which has been constructed to fit the truely 
typical galaxy at each redshift. For example, if color is an indicator of age and/or mctallicity, then the 
results above suggest that our sample contains the oldest and/or most metal rich part of the high redshift 
population. If the objects which satisfy our apparent magnitude limit are in fact, older than the typical high 
redshift galaxy, then it may be that those objects are similar in age to the average object at lower redshifts 
in our sample. If so, then we are better-off using a nonevolving K-correction even though the higher redshift 
sample as a whole is younger. None of the results presented in the main text change drastically if we use 
non-evolving rather than evolving K-corrections. 

To decide which K-corrections to use, we computed the color-magnitude and color-cr relations discussed 
in Paper IV using both the Coleman et al. (1980) template and the Bruzual & Chariot (2003) models. 
The slopes of the mean relations, and the scatter around the mean, remained approximately the same for 
both K-corrections, so we have chosen to not show them here. This suggests that our ignorance of the true 
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Fig. 15. — Apparent g* — r* colors of simulated galaxies in mock catalogs of a passively evolving population; 
the galaxies at higher redshift are younger and, in their rest-frame, bluer. Top panel shows the expected 
distribution of observed colors if there were no magnitude limit; bottom panel shows the effect of imposing 
the same magnitude limit as in our SDSS sample. Solid curves (same in both panels) show the trend of 
observed color with redshift of this evolving population. Although the smooth curve describes the complete 
catalog well, it is substantially bluer than the subset of objects which are included in the magnitude limited 
sample. The difference between the curves and our magnitude limited mock catalogs is similar to that 
between the curves and the data (see previous Figure) , suggesting that the colors of the galaxies in our data 
are evolving similarly to how we assumed in our simulations. 
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K-correction does not strongly compromise our conclusions about how color correlates with magnitude and 
velocity dispersion. Conclusions about evolution, however, do depend on the K-correction. 

The color-a relations constructed using either Kg™ 1 or Kq™J w show evidence for evolution. However, 
K B ™' yields evolution in g* — r* of 0.04 mags, and in r* — i* of 0.07 mags, whereas Kg$ w has changes of 
0.07 and 0.03, respectively. Thus, K B ™' suggests that the evolution in r* — i* is larger than in g* — r* . This 
is not the expected trend; the g* — r* and r* — i* wavelength baselines are about the same, so one expects 
more of the evolution to come in at the bluer color. Using Kg^ w instead suggests that most of the evolution 
is in g* — r* , which is more in line with expectations. 

We also tried K-corrections from Fukugita et al. (1995). At low redshifts, the predicted early-type 
colors are redder than those in our sample, the predicted SO colors are bluer, and the differences depend 
on redshift. A straight average of the two is an improvement, although the resulting low redshift g* — r* is 
red by 0.05 mags. If we shift by this amount to improve the agreement at low redshifts, then the observed 
g* — r* colors at z = 0.25 are redder than the predicted no evolution curves by about 0.2 mags. This is larger 
than the offset we expect for the selection effect introduced by the magnitude limit, so we decided against 
presenting further results from these K-corrections. 

B. Velocity dispersion: methods and measurements 

This Appendix describes how we estimated the line-of-sight velocity dispersions a for the sample of 
galaxies selected for this paper. 

Estimates of a arc limited by the instrumental dispersion and resolution. Recall that the instrumental 
dispersion of the SDSS spectrograph is 69 km s _1 per pixel, and the resolution is ~ 90 km s _1 . In addition, 
the instrumental dispersion may vary from pixel to pixel, and this can affect measurements of a. These 
variations are estimated for each fiber by using arc lamp spectra (upto 16 lines in the range 3800-6170 A 
and 39 lines between 5780-9230 A). An example of the variation in instrumental dispersion for a single fiber 
is shown in Figure 16. The figure shows that a simple linear fit provides a good description of this variation. 
This is true for almost all fibers, and allows us to remove the bias such variations may introduce when 
estimating galaxy velocity dispersions. 

A number of methods for making accurate and objective velocity dispersion measurements as have been 
developed (Sargent et al. 1977; Tonry & Davis 1979; Franx, Illingworth & Hcckman 1989; Bender 1990; Rix 
& White 1992). These methods are all based on a comparison between the spectrum of the galaxy whose 
velocity dispersion is to be determined, and a fiducial spectral template. This can either be the spectrum of 
an appropriate star, with spectral lines unresolved at the spectra resolution being used, or a combination of 
different stellar types, or a high S/N spectrum of a galaxy with known velocity dispersion. In this work, we 
use SDSS spectra of 32 K and G giant stars in M67 as stellar templates. 

Since different methods can give significantly different results, thereby introducing systematic biases 
especially for low S/N spectra, we decided to use three different techniques for measuring the velocity 
dispersion. These are 1) the cross-correlation method (Tonry & Davis 1979); 2) the Fourier-fitting method 
(Tonry & Davis 1979; Franx, Illingworth & Heckman 1989; van der Marel & Franx 1993); and 3) a modified 
version of the direct-fitting method (Burbidge, Burbidge & Fish 1961; Rix & White 1992). Because a galaxy's 
spectrum is that of a mix of stars convolved with the distribution of velocities within the galaxy Fourier 
space is the natural choice to estimate the velocity dispersions — the first two methods make use of this. 
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Fig. 16. — Variation of instrumental dispersion over the range in wavelengths used to measure velocity 
dispersions later in this paper. Solid line shows a linear fit. 

However, there are several advantages to treating the problem entirely in pixel space. In particular, the 
effects of noise are much more easily incorporated in the pixel-space based direct-fitting method. Because 
the S/N of the SDSS spectra are relatively low, we assume that the observed absorption line profiles in 
early-type galaxies are Gaussian (see Rix & White 1992 and Bender, Saglia & Gerhard 1994 for a discussion 
of how to analyze the line profiles of high S/N spectra in the case of asymmetric profiles). 

It is well known that all three methods have their own particular biases, so that numerical simulations 
must be used to calibrate these biases. In our simulations, we chose a template stellar spectrum measured 
at high S/N, broadened it using a Gaussian with rms o~i n put, added Gaussian noise, and compared the input 
velocity dispersion with the measured output value. The first broadening allows us to test how well the 
methods work as a function of velocity dispersion, and the addition of noise allows us to test how well the 
methods work as a function of S/N . 

The best-case scenario is one in which there is no 'template mismatch': the spectrum of the template 
star is exactly like that of the galaxy whose velocity dispersion one wishes to measure. Figure 17 shows the 
fraction of systematic bias associated with each of the different methods in this best-case scenario. Slightly 
more realistic simulations, using a combination of stellar spectra as templates, were also done. The results 
are similar to those shown in Figure 17. With the exception of the cross- correlation method at low (a < 100 
km s _1 ) velocity dispersion, the systematic errors on the velocity dispersion measurements appear to be 
smaller than ~ 3%. 

Although the systematics are small, note that the measured velocity dispersion is more biased at low 
velocity dispersions (a < 100 km s _1 ). For any given S/N and resolution, there is a lower limit on the 
velocity dispersion measurable without introducing significant bias. Since the S/N of the SDSS spectra is 
not very high (see, e.g., Figure 18), and because the instrumental resolution is ~ 90 km s _1 , we chose 70 
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Fig. 17. — Systematic biases in the three methods used to estimate the velocity dispersion. Solid, dashed and 
dotted lines show the biases in the Fourier- fitting, direct-fitting and cross- correlation methods, as a function 
of velocity dispersion and signal-to-noise. 
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Fig. 18. — Comparison of the various methods used to estimate the velocity dispersions; the agreement is 
quite good, with a scatter of about five percent. Most of our spectra have S/N ~ 15, with approximately 
exponential tails on either side of this mean value. 
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Fig. 19. — Distribution of errors as a function of S/N (top) and comparison of estimates from repeated 
observations (bottom). Both panels suggest that, when the S/N > 15 then the typical error on an estimated 
velocity dispersion is S\og 10 a < 0.04. 
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km s as a lower limit. Figure 18 shows a comparison of the velocity dispersion estimates obtained from 
the three different methods for the galaxies in our sample. The median offsets are not statistically significant 
and the rms scatter is ~ 0.05. On the other hand, the top panels suggest that the cross-correlation method 
sometimes underestimates the velocity dispersion, particularly at low S/N. 

We evaluate the dependence of the velocity dispersion on the wavelength range by fitting the spectra in 
different intervals: 4000 — 5800 A which is the usual wavelength range used in the literature; 3900 — 5800 A 
to test the effect of including the Ca H and K absorption lines (e.g., Kormendy 1982); 4000 — 6000 A to test 
the effect of including the NaD line (e.g., Dressier 1984); and 4000 - 7000 A and 4000 - 9000 A to test the 
effect of including longer wavelengths. The velocity dispersion obtained with the Ca H and K is ~ 2% larger 
than that obtained using the standard wavelength region 4000 — 5800 A, and the rms difference between the 
three different methods increases to <~ 7%. Including the NaD line increases the velocity dispersion by <~ 3% 
but does not increase the scatter between the different methods. Using the wavelength range 4000 - 7000 A 
only provides velocity dispersions which are ~ 3% larger than the values obtained if only 4000 — 5800 A 
range is used. On the other hand, in this wavelength region, the different methods (and measurements from 
repeated observations) are in better agreement; the scatter is ~ 8% smaller than in the 4000 — 5800 A region. 
In the range 4000 — 9000 A, the velocity dispersion estimates increase by ~ 7%. This last effect is probably 
due to the presence of molecular bands in the spectra of early-type galaxies at long wavelengths (i.e., to the 
presence of cool stars). Furthermore, the scatter in this wavelength region increases dramatically (<~ 15%). 
Presumably this is due to the presence of higher sky-line residuals and lower S/N . 

The estimated velocity dispersion we use in the main text are obtained by fitting the wavelength range 
4000 — 7000 A and then using the average of the estimates provided by the Fourier-fitting and direct-fitting 
methods to define what we call a cst - We do not use the cross-correlation estimate because of its behavior at 
low S/N as discussed earlier. 

The top panel of Figure 19 shows the distribution of the errors on the velocity dispersion as a function 
of the S/N of the spectra. The errors for each method were computed by adding in quadrature the statistical 
error due to the noise properties of the spectrum, and the systematic error associated with the template 
and galaxy mismatches. The final error on a cs t is got by adding in quadrature the errors on the two 
estimates (i.e., the Fourier- fitting and direct-fitting) which we average. The resulting errors range from 
0.02 < 5 log 10 a < 0.06 dex, depending on the S/N of the spectra, with a median value of 0.03 dex. 

A few galaxies in our sample have been observed more than once. The bottom panel shows a comparison 
of the velocity dispersion estimates from multiple observations. The scatter between different measurements 
is ~ 0.04 dex, consistent with the amplitude of the errors on the measurements. 

C. Velocity dispersion: profiles and aperture corrections 

The SDSS spectra measure the light within a fixed aperture of radius 1.5 arcsec. Therefore, the estimated 
velocity dispersions of more distant galaxies arc affected by the motions of stars at larger physical radii than 
for similar galaxies which are nearby. If the velocity dispersions of early- type galaxies decrease with radius, 
then the estimated velocity dispersions (using a fixed aperture) of more distant galaxies will be systematically 
smaller than those of similar galaxies nearby. 

We have not measured the velocity dispersion profiles cr(r) of any of the galaxies in our sample, so we 
cannot correctly account for this effect. If we assume that the galaxies in our sample are similar to those for 
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Fig. 20. — Velocity dispersions of galaxies as a function of redshift. Top panel shows the estimated velocity 
dispersion, and bottom panel shows the values after correcting the estimate as described in the main text 
(equation 1). Different symbols show the result of averaging over volume-limited subsamples (same as in 
Figure 12) of galaxies having approximately the same luminosities and effective radii at each redshift. (Error 
bars show the rms scatter around this mean value.) The mean trends with redshift can be used to infer how, 
on average, the velocity dispersion changes with distance from the centre of the galaxy, and how this change 
depends on luminosity and effective radius. 
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Fig. 21. — As for previous figure, but for galaxies with larger radii. 
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which velocity dispersion profiles have been measured, then we can use the published a(r) curves to correct 
for this effect. This is what equation (1) in Section 3.3.1 does. 

An alternative procedure can be followed if evolution effects are not important for the velocity dispersions 
in our sample. To illustrate the procedure, the galaxies in each of the volume-limited subsamples shown 
in Figure 12 were further classified into small bins in effective physical radius (i.e., R a in kpc/h, not r Q 
in arcsec). Figures 20-21 show the result of plotting the velocity dispersions of these galaxies versus their 
redshifts. Since the galaxies at low and high rcdshift are supposed to be similar, any trend with rcdshift 
can be used to infer and average velocity dispersion profile, and how the shape of this profile depends on 
luminosity and effective radius. In this way, the SDSS data themselves can, in principle, be used to correct 
for the effects of the fixed aperture of the SDSS spectrograph. 

In practice, because there is substantial scatter in the velocity dispersions at fixed luminosity and size 
(Paper II shows that this scatter is about 14%), the trends in the present data set are relatively noisy. When 
the dataset is larger, it will be worth returning to this issue. For now, because the corrections are small 
anyway, we have chosen to use equation (1) to correct the velocity dispersions. Nevertheless, curves like 
those presented above provide a novel way to study the velocity dispersion profiles of early-type galaxies. 



Let £ denote the covariance matrix of the errors in our estimates of the the absolute magnitude M, the 
half light radius R = log 10 R 0l and the velocity dispersion V = log 10 a: 



The elements of the error matrix £ are obtained as follows. 

The photometric pipeline estimates the size rd ov and the apparent magnitude rridev from the same fitting 
procedure. As a result, errors in these two quantities are correlated. Let e r denote the error in log 10 rd cv , and 
e m the error in TOdcv The correlation means that we need three numbers to describe the errors associated 
with the fitting procedure, (e r e r ), (e m e m ), and (e m e r ), but the pipeline only provides two. The error output 
by the pipeline in rd cv , is correctly marginalized over the uncertainty in TOdcv, so it is essentially (e r e r ). On 
the other hand, the quoted error in TOdcv, say (ep hoto ) is really (e m e m ) — (e m e r ) 2 /(e r e r ). To estimate the 
values of (e m e r ) and (e m e m ) which we need, we must make an assumption about the correlation between 
the errors. 

Fortunately, this can be derived from the fact that, for a wide variety of galaxy profile shapes, the 
quantity £ = e r — ae^, with a « 0.3, has a very small scatter (e.g. Saglia ct al. 1997). Here ji = 
TOdcv + 51og 10 rd C v + 2.5 log 10 (27r) is the surface brightness, and e M is the error in the surface brightness. As 
a result, 
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(Saglia et al. 1997). This means that (e r £) = (£ 2 )/(l + a 2 ), so that 
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The final equality shows that the error output from the pipeline provides an estimate of (£ 2 ) which we can 
insert into our expressions for (e m e m ), (e m e r ), and (e^e^). (Notice that if (£ 2 ) -C (e r e r ), then it would be 
a good approximation to set (e 2 hoto ) ~ (£, 2 )/a 2 . Since this is not always the case for our dataset, we must 
solve the quadratic.) Once this has been done, we set 
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That is, we compute the error in the absolute magnitude by assuming that there are no errors in the 
determination of the redshift (and K correction!) which would otherwise propagate through. 

Subsequent papers will focus almost exclusively on the circularly averaged radius R defined in Sec- 
tion 3.2.1. The errors on it are given by adding the errors in the size r<i cv to those which come from the 
error on the shape b/a. We assume that the errors in b/a are neither correlated with those in log 10 r^cv nor 
with those in the absolute magnitude. Finally, we assume that errors in magnitudes are not correlated with 
those in velocity dispersion, so (ey M ) is set to zero, and that errors in size and velocity dispersion are only 
weakly correlated because of the aperture correction we apply. Here {e v e v ) is the error in what was called 
log 10 (Test in the main text. 

These error estimates are presented in Tables 4 and 4, where we have set Sm p hoto — \/( e m)' $r — \/( e r)' 
and SM — £mm- 



